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ABSTRACT 

In this paper, we discuss the dynamical issues of quantum computation. We demon- 
strate that fast wave function oscillations can affect the performance of Shor's quantum 
algorithm by destroying required quantum interference. We also show that this destruc- 
tive effect can be routinely avoided by using resonant-pulse techniques. We discuss the 
dynamics of resonant pulse implementations of quantum logic gates in Ising spin sys- 
tems. We also discuss the influence of non-resonant excitations. We calculate the range 
of parameters where undesirable non-resonant effects can be minimized. Finally, we 
describe the "27r/c-method" which avoids the detrimental deflection of non-resonant 
qubits. 
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I. Shor's quantum algorithm — the simplest example 

Successful design of quantum computers requires a thorough understanding of the 
time evolution of quantum qubits. The current interest in quantum computation was 
stimulated by Shor [[[]] who invented the quantum algorithm for prime factorization of 
integers. Shor's algorithm has three main steps. (See (7) through (11) below.) We 
discuss here the simplest example of Shor's algorithm and use it in the next section 
to demonstrate the role of fast oscillations, e Suppose we want to factor the smallest 
composite number, 4. According to Shor's method, we choose the only coprime number, 
3. A quantum computer will then compute the periodic function, 

y(x) = 3 X (mod 4), (1) 

which is the remainder after division of 3 X by 4. This function is, 

y(0) = l, 2/(1) = 3, y(2) = 1, y(3) = 3,... (2) 

Then, the quantum computer must find the period of this function, T = 2, in our case. 
The factor of 4 can be found as the greatest common divisor (GCD) of two numbers: 
(z + 1) and 4 or (z — 1) and 4, where z = 3 T//2 . In our case, 

GCD(z -1,4) = GCD(2, 4) = 2, (3) 

gives the desired factor of 4. 

The simplest quantum computer which computes the function (1) and finds its 
period, T, has 4 quantum bits (qubits): the two left-most qubits represent the number 
"x" , and the two right-most qubits represent "y" , 

x = 2m 1 + m , y = 2n\ + n , m iy ni = 0,1, (i = 0, 1). (4) 

As an example, the values x — 1 and y = 3 will be represented as, 

\mimo, niTio) = |01, 11). (5) 

We use the Dirac notation for the quantum states: |0) for the ground state and |1) for 
the excited state. 

Assume that initially all four qubits are in their ground states. The wave function 
of the system is, 

^o=|00,00). (6) 

Following Shor's idea, the quantum computer must carry out three main unitary trans- 
formations: 

1) It creates a superposition of all possible values of "x", 

® 1 = (|00, 00) + |01, 00) + |10, 00) + |11, 00)). (7) 
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2) It computes the function, y(x), in (1) using a digital algorithm, for all values of "x" 
simultaneously, 

# 2 = ~(|00,01> + |01,11) + |10,01> + |11,11>). (8) 

3) It performs a discrete Fourier transform (DFT) for each value of "x" , 

\x) -> -j2e 2 ™ kx/4 \k). (9) 

^ fc=0 

For example, 

|00,01) -> -(|00,01) + |01,01) + |10,01) + |11, 01)). (10) 
2 

As a result, 

^ 3 = ^(|oo, 01) + |00, 11) + |10, 01) - |10, 11)). (11) 

Measuring the value of "x" one gets either: x = X\ = or x = X2 = 2. The ratio D/x2 
(where D = 4 is the number of all possible values of "x" ) is the period of the periodic 
function y(x): T = Djx-i = 2 @. 

II. Fast oscillations 



Note that Shor's algorithm is described above as a sequence of instantaneous unitary 
transformations. What happens if each unitary transformation takes a finite time? In 
time, t, in quantum mechanics, each stationary state, \k), of the quantum system 
acquires its own phase, —Ekt/h. This phase is usually responsible for fast oscillations 
of the wave function. The slow dynamics is caused by small perturbations which are 
responsible for transitions between the energy levels. Generally it is not obvious that 
the DFT will provide the desired constructive and destructive interferences if one takes 
into account the fast oscillations of wave function during the finite time intervals of 
quantum transformations and possible time delays between these transformations. 

Suppose, for simplicity, that there are finite time delays between three main unitary 
transformations described in section I, while the time duration of each transformation 
is, as before, infinitely small. Then, immediately before the second transformation, 
each state in the superposition (7) will possess the "natural" phase factor, 

q? 1 = i|j00,00)e-^ OOTl + |01, 00> e -^ 10 " 1 + [10, 00)e~^ aon + 1 11, OO^"^ 30 ^ , (12) 

where we put h — 1; E a b is the energy corresponding to the state: x — a, y — b; T\ 
is the time delay between the first two transformations. The second transformation 
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generates new states, and the "natural" connection between the state and the phases 
vanishes. Just before the third transformation we have, 

^ 2 = {|00, 01) expH^ooTi - iE 01 r 2 ) + |01, 11) exp(-^ 10 r 1 - iE 13 r 2 )+ (13) 

1 10, 01) exp(-iE 20 r 1 - iE 21 t 2 ) + |11, 11) exp(-iE 30 ri -iE 33 T 2 )}, 

where r 2 is the time delay between the second and the third unitary transformations. 
It is easy to check that the desired interference does not occur after the DFT. For 
example, the state (1/2)|00,01) in the superposition (11) is described now by two 
terms, 

-|00,01)[exp(-i£ on - iE 01 T 2 ) + exp(-iE 20 Ti - iE 21 r 2 )}. (14) 
The state |01, 01) (which vanishes if Ti = r 2 = 0) now survives, and has the form, 

i|01,01)[exp(-iE 00 Ti - iE 01 r 2 ) - exp(-iE 20 ri - iE 21 r 2 )]. (15) 

The reason for breaking Shor's algorithm in this example is that the "like" terms 
corresponding to the same state, carry the "history" of their origination. For example, 
the term, 

-\00, 01) eM-iEoon - iE 01 r 2 ), (16) 
in (14) originated from the transformations, 

|00, 00) -> |00, 00) ^ |00, 01) A |00, 01). (17) 

The term, 

-|00,01)exp(-^ 20 r 1 -iE 21 r 2 ), (18) 
originated from the transformation, 

|00, 00) -> 1 10, 00) ^ |10, 01) A |00, 01). (19) 

One can see that the "phase memory" includes the energies connected with the process 
of formation of the corresponding terms. The question arises: "How can one carry out 
a quantum computation with finite-time intervals between pulses and not destroy the 
desired interferences?" 

The natural way to solve this problem is to "produce" the "natural" phase, —E ab At, 
when the term corresponding to x = a and y = b is generated in the process of quantum 
computation. (Here, At, is the time interval between the beginning of the first trans- 
formation and the end of the transformation which generates this state, \x — a, y — b).) 
We now discuss how this can be done routinely for a resonant technique such as NMR. 
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III. Resonant pulse implementation of quantum unitary transformations 



Recently, two basic approaches have been used for implementing unitary trans- 
formations for quantum computation. The most popular approach originated from 
the paper by Lloyd who proposed using resonant electromagnetic pulses to perform 
any required transformation ||. These pulses drive the quantum transitions between 
the corresponding energy levels in a system of weakly interacting particles. Following 
Lloyd's idea, we proposed in || using Ising spin systems for quantum computation. 
For such systems, the "digital states" corresponding to specific values of x and y, are 
true eigenstates of the Hamiltonian including interaction between spins. This idea has 
been implemented in nuclear spin systems in liquids, which are closely approximated 
by Ising systems ^ Ej. It was found in || [J that quantum computation can be carried 
out on a statistical ensemble of such systems at room temperature. There are also a 
few other important "resonant" proposals including a well-known idea of the ion trap 
quantum computer 0. 

The second basic approach does not require resonant pulses. The earliest proposals 
of this kind were reviewed in ||, and the most recent idea was suggested in ||. Here, 
we are considering only the "resonant" proposals. In this section, we discuss how the 
"natural" phase, —E a bAt, can be routinely generated in a quantum computer using 
resonant pulses. 

We assume that there is a required sequence of resonant pulses to produce a quan- 
tum computation. Each pulse induces a selective transition between two definite states 
of the system. We denote the state of the whole system by a single index, e.g. \n). The 
transition between two states, \k) and \n), can be described by the following effective 
Hamiltonian: 



H = E k \k)(k\+E n \n)(n\-^ 



A»*v)\k)(n\ + e~ i{ult+<p) \n){k\ 



(20) 



Here E k and E n are the energies of the corresponding states; Q is the Rabi frequency; 
us = u n k = E n — Ek is the resonant frequency, and <p is the phase of the resonant pulse. 
We assume that one imposes the following resonant fields, 



h + = h exp(-iu nk t) , 



(21) 



for each resonant transition required for the quantum computation. Suppose that the 
j-th pulse induces a transition from state \k) to state |n). We assume that at the 
beginning of the pulse, the amplitude, C k , of the state \k) has the "natural" phase 
factor, 

exp(-iE k tj-i), tj-x = n + T 2 + ... + Tj-i, (22) 

where r k is the duration of the fc-th pulse. We also assume that there is no delay between 
two consecutive pulses. (This last assumption is not crucial for our consideration.) The 
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Schrodinger equation for the amplitudes C k and C n has the form, 

rC^EA-h^C,, i6.-E.C.-\ a r^C k . (23) 
At the end of the j-th pulse (t = tj = tj-i + Tj), we have, 

C k (tj) = exp[i(E k tj-i - E k tj)} coaaijCk{tj-i), ctj = Qtj/2, (24) 

C n {tj) = exp[i(7r/2 - <p + E k t^i - E n tj)\ saxajCk{tj-i). 
Using (22), we can rewrite these expressions in the form, 

Ck(tj) = \Ck(tj-i)\ exp(-iE k tj) cosaj, (25) 

C n (tj) = |C fc (tj_i)| sino! J exp[(i(7r/2 - <p)) exp(-iE n tj). 

One can see from (25) that the new state, |n), has automatically acquired the "natural" 
phase, —E n tj. The factor, exp[z(7r/ 2 — </?)], in (25) describes the "standard" phase shift 
which can be eliminated by choosing: ip = tt/2. 



IV. Non-resonant action of resonant pulses 



A major obstacle for the resonant pulse implementation of quantum computation is the 
non-resonant effect of resonant pulses on some qubits. This action can be completely 
avoided only if each qubit is isolated from all other qubits. This is not possible or 
desirable in reality. One can minimize non-resonant effects by providing a significant 
difference between the resonant frequencies. We shall now discuss this possibility us- 
ing the example of a two-qubit quantum Control-Not gate which is described by the 
operator, 

CN = |00)(00| + |01)(01| + 1 10> <1 1 1 + |11)(10|. (26) 

This gate changes the value of the right (target) qubit if the left (control) qubit is in 
the excited state. The importance of the CN gate became evident when it was shown 
that any unitary transformation can be implemented using combinations of one-qubit 
rotations and two-qubit CN gates ||10||. A single 7r-pulse implementation of the CN gate 



was first suggested for coupled quantum dots [p| . We have studied this implementation 
using Ising spin systems [f[2]j-[rT|. We investigated the dynamics of a single 7r-pulse 



CN gate in the reference frame rotating with frequency, lu, of the circularly polarized 
magnetic field, which induces a resonant transition of a target spin, if the control spin 
is in the excited state. The effective Hamiltonian in the rotating frame has the form 

[0: 



n 



E 



2 J kn I k I^ 



(27) 



6 



where Uk is the resonant frequency of a single fc-th spin and J^ n are the Ising interaction 
constants. The utility of this Hamiltonian for liquid NMR experiments was tested and 
confirmed in [pl|. We have shown |14|, [17) that a single pulse performs the CN gate 
operation, 

CM = |00)(00| + |01)(01| +z|10)(ll| + i|ll)(10|, (28) 

accompanied by a 7r/2 phase shift for the transformation of the target qubit. Numerical 
calculations for a well-separated frequencies, 



u x = 500, to 2 = 100, J = 5, 



uj 2 -J = 95, fii = 0.5, fi 2 = 0.1, (29) 



confirmed that non-resonant effects are negligible for the single-pulse implementation 
of the CN gate [[T^, [T7[] . For example, the initial superpositional state, 



^o = v^3|00) + VK2\01) + -J=|10) + -±=|11>, 

transforms as a result of the CM gate (28) to the state (See Fig 25.2 in |l 
m x = y/03\00) + VM\0l) + (i'/V6)|10) + (i/v^lll). 



(30) 



(31) 



We have also investigated the dynamics of the CN gate for an ensemble of four-spin 
molecules at room temperature: k B T » TiUk Hl5l-[pJ|. Following the idea suggested in 
[§], we implemented quantum transformations for four active states, \00ij), = 0, 1) 
which can be prepared initially in the ground state, |0000). The corresponding density 
matrix has the form, 

p(t) = E/lQ + p A (t), (32) 
where pa is the following deviation density matrix: 



32k B T 



E r n , P (t)\n){p\ + J2 b nAt)\n){p\ 

n,p=0 



(33) 



The first sum in (33) describes the dynamics of the 16 density matrix elements, r ntP (t), 
corresponding to the dynamics of a superpositional "active" state, J2ij Cij\00ij). The 
second sum describes the dynamics of all the remaining 240 density matrix elements. 
These matrix elements should not change significantly during the dynamical process. 
Here we use a "single-index" decimal notation for each four-spin state, 



\n) = \psij), n = j + 2i + 2 2 s + 2 3 p. 



(34) 



The initial conditions for the matrix elements b njP are: 65,5 



h,6 = b 7 ,7 = b s ,s = 1/2, 

64,4 = &9,9 = &10,10 — &11.11 = — l/2, 6l2,12 = — 1) &13.13 = &14,14 = &15.15 = 0, b n)P = for 

n ytz p. Our numerical calculations confirmed that single-pulse implementation of the 
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CN gate performs well if the frequencies, uuk, are well separated and the values of the 
Ising constants are not too small. For example, for the values of parameters, 



UJ k = UJ + 100 X k, Jkn = J = 10, fife = fl = 0.1, 



(35) 



we have demonstrated the "complementary" CN gate which drives the target spin if 
the control spin is in the ground state. For the initial superpositional state (30), after 
the action of the complementary CN gate, one creates the following wave function for 
the effective pure two-spin system of two "active" qubits, 



= iv / cT2|00> + zVO3|01) + 4=|!0) + ^l 11 )- 

V3 v6 



(36) 



It is easy to check that the reduced density matrix, r, corresponding to the initial state 
(30) has the form: 



r 



^0.3000 0.2449 0.3162 0.2236^ 

0.2449 0.2000 0.2582 0.1826 

0.3162 0.2582 0.3333 0.2357 

V 0.2236 0.1826 0.2357 0.1666/ 



(37) 



The density matrix corresponding to the wave function (36) is, 



n 



( 0.2000 0.2449 z0.2582 i0.1826^ 

0.2449 0.3000 i0.3162 i0.2236 

-i0.2582 -z0.3162 0.3333 0.2357 

\-i0.1826 -i0.2236 0.2357 0.1666 / 



(3* 



The numerical solution of the equations of motion corresponding to the density matrix, 
Pa(£), in (33) showed less than 0.5% deviations of r„ iP from the values in (38). With the 
same accuracy, the coefficients b njP (t) in (33) do not change in time under the action 
of the 7r-pulse. 

Next, we have studied numerically the range of parameters in which non-resonant 
effects remain small (with about 1% deviations from the dynamics corresponding to 
the CN gate) [[Uj. According to our calculations, non-resonant effects do not destroy 
the single-pulse CN gate when Auj/Q > 300, where Au; = u>k+i — u>k- This is valid 
for a wide range of the values of the Ising constant of interaction: 3 < J/Q < 100. 
An analysis of a single-pulse CN gate was also performed in |18| where it was called 
the "Pound- Overhauser implementation of the CN gate" after the Pound-Overhauser 
double resonance in the NMR spectroscopy (on-transition excitation). 

To conclude this section, we mention related topics associated with a CN gate in 
Ising spin systems. The influence of small deviations from resonance and effects of 

121 . Two-pulse implementations of the CN gate have 



"noise" have been studied in 
been considered in QT9L [2(|. In this case, the first short pulse induces a n/2 rotation 



8 



of the target spin around the x-axis of the rotating reference frame. The second pulse 
induces a similar rotation around the y-axis. The delay time between the two pulses 
is n/2J. The action of the second pulse depends on the state of the control spin. 



For example, in experiments []!£], the control proton spin has the resonant frequency 
500MHz; the target carbon spin has the resonant frequency 125MHz; and the Ising 
constant is: J/2tt = 108Hz. The two 7r/2-pulse sequence is more complicated than 
a single 7r-pulse implementation of the CN gate but it takes less time (approximately 
7r/2J). Another two-pulse implementation of the CN gate was called the "Pound- 
Overhauser on-resonance implementation of the CN gate" JT8|- It consists of a long 
x-pulse (with a duration of 7r/(\/2J), and a Rabi frequency of Q — J) which is tuned 
to the frequency between the peaks of the NMR doublet for the target spin, and a 
short 7r/2-pulse (y-pulse) covering the doublet. 

V. Elimination of non-resonant excitation 



The question arises about the opportunities to eliminate non-resonant effects using 
special sequences of resonant pulses. One well-known approach to remove "parts" 
of the Hamiltonian in nuclear spin systems is "selective averaging". This is one of 



the basic approaches for high resolution NMR in solids. (See, for example, [2TJ and 
references therein.) It was suggested recently to use this approach for analog quantum 
computations |E£9. Indeed, using averaging (a proper sequence of pulses) one can 
transform with some accuracy the initial Hamiltonian into the desired Hamiltonian to 
study the evolution of the desired quantum system. In [E2] this was done for two qubits, 



but the method can probably be extended for a larger number of qubits. Problems 
which one faces using this method have been discussed already in the first papers 
H,|2J. (See also JTLj.) 

To the best of our knowledge, no one has proposed eliminating non-resonant part 
of the Hamiltonian using averaging. Instead of averaging, we suggested in [|I^, 17 1 a 
different method (the "27rfc-method"): a 7r- or n/n (n = 2,3,..)- pulse which drives 
the resonant spin must simultaneously be 27rfc-pulse for a non-resonant spin. In this 
case, the non-resonant spin will return to its initial unperturbed state at the end of 
the excited pulse. 

Consider, for example, two spins with a small frequency difference, Aw. If the 
exciting pulse, e.g. a 7r-pulse, is tuned to the frequency of the first spin, we have, 

tir = it, (39) 

where r is the duration of the pulse and Q is the Rabi frequency. For the non-resonant 
spin, the effective field in the rotating frame in frequency units is, 



Vtt 2 + Aw 2 , (40) 
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where we have assumed equal Rabi frequencies for both spins. If, 

u e T = 2nk, k = l,2,..., (41) 

the non-resonant spin is not deflected by the pulse. From (39)-(41) we easily derive the 
required conditions for the Rabi frequency and for the duration of the 7r-pulse which 
allow one to eliminate undesirable non-resonant excitations |R|. (See also Chapter 22 
in 0.), 

a =vwh- T ^ /a - (42) 

(We do not discuss here the phase effect for a 27r-pulse: for a pure quantum state, a 
27r-pulse changes the phase of the state by 7r and a 47r-pulse returns the system into the 
initial quantum state. See, for example, [17].) By manipulating Au, one can eliminate 
non-resonant effects for two or more weakly interacting spins. Assume, for example, 
that the frequencies of spins are, 

u , uji = u +8Q, u 2 = U0 + I6Q, u 3 = co> +24f2, ...,u n = u +8nQ, (n = 1,2,...). (43) 

If one applies a resonant 7r-pulse with any of these frequencies, e.g., u = u 2 , the 
frequency difference, |Au;|, takes the values, 

8Q, 16Q, 24vr, ...,8nQ, (n = l,2, ...). (44) 

According to (40), (41), (44) the corresponding angles of rotation for the non-resonant 
spins are, 



Vfl 2 + AuJ 2 X J = 7T, 



/AuA 2 _ tt\Auj 



+ lm ' 1 = 8?r, 16vr, 24vr, 8nvr, n = l,2, ... . 



(45) 

It is clear from (45) that at the end of the 7r-pulse, all non-resonant spins will return 
to states which are very close to the initial states. 

For Q = 10 3 s _1 , the frequency difference is Au = 8Q which corresponds to a change 
of the external magnetic field, B, for the neighboring protons by: AB s;2x 10~ 5 T. 
Choosing the distance between neighboring protons, Ax = lnm, we obtain an estimate 
for the required magnetic field gradient, 

AB 

«3x 10 4 T/m. (46) 

Ax 

This gradient is easily achievable, for example, in recent Magnetic Resonance Force 
Microscopy experiments [f2~5|| . 

For a 7r/n-pulse (n = 2, 3, ...), instead of (42) we have the generalized condition, 

= —n |A ^' » (A: =1,2,...). (47) 
\2nk) 2 - 1 
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It was suggested in |L8j that the same 27r/c-method be used to provide an "exact" 
single-pulse CN gate in the two-spin system. If a control spin is in the ground state, 
the frequency of the transition for the target spin is ujq + J, where ujq is the target spin 
frequency ignoring Ising interactions. If the control qubit is in the excited state, the 
frequency of the transition is ujq — J. So, in this case, the frequency difference in (42) 
is: Au = 2J. A pulse with frequency uo = ujq — J and with the Rabi frequency and 
duration from (42) represents a resonant 7r-pulse if the control spin is in the excited 
state, and a 27r£;-pulse if the control spin is the the ground state. In the second case, 
the target spin returns to its initial state. Choosing k — 1, one can reduce the duration 
of a single-pulse CN gate to: ^/3ti/2J. 



Conclusions 



In this paper we discuss problems related to the dynamics of quantum computation. 
We presented the simplest demonstration of the dynamics of Shor's algorithm using 
only four qubits. In this example, we showed that fast oscillations can destroy the 
desired quantum interference. We also demonstrated how this effect can be routinely 
eliminated using resonant pulses. We also considered influence of non-resonant effects 
on quantum computation. We have shown that for well-separated frequencies, a non- 
resonant excitation does not prevent a single-pulse implementation of logic gates for 
both pure quantum states and for an ensemble of quantum systems. We have found 
the range of parameters where non-resonant effects remain small. Finally, we have de- 
scribed a u 2irk -method" which allows one to drive resonant spins without deflecting of 
non-resonant spins. Future studies of dynamical problems related to quantum compu- 
tation will be directed to many-qubit systems and the influence of noise in real systems. 



Acknowledgments 



This work was supported by the Department of Energy under contract W-7405-ENG- 
36, and by the National Security Agency. 



11 



References 

[1] P. Shor, Proc. of the 35th Annual Symposium on the Foundations of Computer 
Science, IEEE, Computer Society Press, New York, (1994), p. 124. 

[2] A. Ekert, R. Jozsa, Rev. Mod. Phys., 68, (1996) 733. 

[3] S. Lloyd, Science, 261, (1993) 1569. 

[4] G.P. Berman, G.D. Doolen, G.D. Holm, V.I. Tsifrinovich, Phys. Lett. A, 193 
(1994) 444. 

[5] N.A. Gershenfeld, I.L. Chuang, Science, 275, (1997) 350. 

[6] D.G. Cory, A.F. Fahmy, T.F. Havel, Proc. Natl. Acad. Sci. USA, 94 (1997) 1634. 

[7] J.I. Cirac, P. Zoller, Phys. Rev. Lett, 74, (1995) 4091. 

[8] R. Landauer, Philos. Trans. R. Soc. London, 353, (1995) 367. 

[9] D. Loss, D.P DiVincenzo, Phys. Rev. A, 57, (1998) 120. 

[10] A. Barenco, C.H. Bennett, R. Cleve, D.P. DiVencenzo, N. Margolus, P. Shor, T. 
Sleator, J.A. Smolin, H. Weinfurter, Phy. Rev. A, 52, (1995) 3457. 

[11] A. Barenco, D. Deutsch, A. Ekert, R. Jozsa, Phys. Rev. Lett, 74, (1995) 4083. 

[12] G.P. Berman, D.K. Campbell, G.D. Doolen, G.V. Lopez, V.I. Tsifrinovich, Physica 
B, 240 (1997) 61. 

[13] G.P. Berman, D.K. Campbell, V.I. Tsifrinovich, Phys. Rev. B, 55, (1997) 5929. 

[14] G.P. Berman, G.D. Doolen, G.V. Lopez, V.I. Tsifrinovich, |quant-ph/9802013| . 

[15] G.P. Berman, G.D. Doolen, G.V. Lopez, V.I. Tsifrinovich, Phys. Rev. B, 58, 
(1998) 11570. 

[16] G.P. Berman, A.R. Bishop, G.D. Doolen, G.V. Lopez, V.I. Tsifrinovich, LA-UR- 
98-4632, submitted to Phys. Rev. B, 1999. 

[17] G.P. Berman, G.D. Doolen, R. Mainieri, V.I. Tsifrinovich, Introduction to Quan- 
tum Computers, World Scientific Publishing Company, 1998. 

[18] D.G. Cory, A.E. Dunlop, T.F. Havel, S.S. Somaroo, W. Zhang, The effective 
Hamiltonian of the Pound- Overhauser controlled-NOT gate, submitted for publi- 
cation, 1999. 



12 



[19] I.L. Chuang, N.A. Gershenfeld, M.G. Kubinec, D.W. Leung, Proc. Royal Soc. 
London Ser. A- Math. Phys. Engin. Sciences, 454, (1998) 447. 

[20] D.G. Cory, M.D. Price, T.F. Havel, Physica D, 120, (1998) 82. 

[21] U. Haeberlen, High Resolution NMR in Solids, Academic Press, 1976. 

[22] S. Somaroo, C.H. Tseng, T. Havel, R. Laflamme, D.G. Cory, Quantum simulations 
on a quantum computer, submitted for publication, 1999. 

[23] J.S. Waugh, C.H. Wang, Phys. Rev., 162, (1967) 209. 

[24] J.S. Waugh, C.H. Wang, L.M. Huber, R.L. Void, J. Chem. Phys., 48, (1968) 662. 

[25] P.C. Hammel, Z. Zhang, G.J. Moore, M.L. Roukes, J. Low Temp. Phys., 101, 
(1995) 59. 



13 



